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Abstract. Boundary value problems for integrable nonlinear evolution PDEs for- 
mulated on the finite interval can be analyzed by the unified method introduced 
by one of the authors and used extensively in the literature. The implementation 
of this general method to this particular class of problems yields the solution in 
terms of the unique solution of a matrix Riemann-Hilbert problem formulated in 
the complex fc-plane (the Fourier plane), which has a jump matrix with explicit 
(x, t)-dependence involving six scalar functions of k, called spectral functions. Two 
of these functions depend on the initial data, whereas the other four depend on all 
boundary values. The most difficult step of the new method is the characterization 
of the latter four spectral functions in terms of the given initial and boundary data, 
i.e. the elimination of the unknown boundary values. Here, we present an effective 
characterization of the spectral functions in terms of the given initial and boundary 
data. We present two different characterizations of this problem. One is based on 
the analysis of the so-called global relation, on the analysis of the equations ob- 
tained from the global relation via certain transformations leaving the dispersion 
relation of the associated linearized PDE invariant, and on the computation of the 
large k asymptotics of the eigenfunctions defining the relevant spectral functions. 
The other is based on the analysis of the global relation and on the introduction 
of the so-called Gelfand-Levitan-Marchenko representations of the eigenfunctions 
defining the relevant spectral functions. We also show that these two different char- 
acterizations are equivalent and that in the limit when the length of the interval 
tends to infinity, the relevant formulas reduce to the analogous formulas obtained 
recently for the case of boundary value problems formulated on the half-line. 
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1. Introduction 

This is the third in a series of papers addressing the most difficult problem in the 
analysis of integrable nonlinear evolution PDEs, namely the problem of expressing 
the so-called spectral functions in terms of the given initial and boundary conditions. 
In [9] this problem was analyzed for the case of the half-line. In |10j the same problem 
was also analyzed with the aid of the so-called Gelfand-Levitan-Marchenko (GLM) 
representations. Here we analyze this problem for the case of the finite interval. 

We refer the interested reader to [9] for an introduction to the unified method 
of [H El Ej and for a discussion of the difference between linearizable versus non- 
linearizable boundary value problems. Here, we only note that the unified method 
expresses the solution q{x,t) of an integrable evolution PDE in terms of an integral 
formulated in the complex fc-plane. This representation is similar to the integral 
representation obtained by the new method for the linearized version of the given 
nonlinear PDE, but it also contains the entries of a certain matrix-valued function 
M(x,t,k), which is the solution of a matrix Riemann-Hilbert (RH) problem. The 
main advantage of the new method is the fact that this RH problem involves a jump 
matrix with explicit (x, t)-dependence, uniquely defined in terms of the spectral func- 
tions. For the problem on the interval, there exist six spectral functions denoted by 
{a{k),b{k), A{k), B{k), A{k), B{k)}. The functions {a{k),b{k)} are defined in terms 
of the initial data qoix) = q{x,0) via a system of linear Volterra integral equations; 
the functions {A{k),B{k)} and {A{k),B{k)} are defined in terms of the boundary 
values at X = and x = L respectively, also via systems of linear Volterra integral 
equations. However, the integral equations defining {A{k),B{k)} and {A{k), 13{k)} 
involve all boundary values, whereas for a well-posed problem only a subset of these 
boundary values can be prescribed as boundary conditions. Thus, the complete so- 
lution of a concrete initial-boundary value problem requires the characterization of 
{A(k), B(k), A{k), B{k)} in terms of the given initial and boundary conditions. For 
example, for the Dirichlet problem of the NLS on the interval, < x < L, it is nec- 
essary to characterize {A{k), B{k),A{k),B{k)} in terms oiqo{x), go{t) := q{0,t), and 
hoit)=qiL,t). 

A characterization of the spectral functions is called effective if it fulfills the fol- 
lowing requirements: (a) In the linear limit, it yields an effective solution of the 
linearized boundary value problem, i.e. it yields a solution in the form of an integral 
which involves the transforms of the given initial and boundary conditions. (6) For 
'small' boundary conditions, it yields an effective perturbative scheme, i.e. it yields an 
expression in which each term can be computed uniquely in a well-defined recursive 
scheme. 

The effective characterization presented here is based on the construction of the 
generalized Dirichlet to Neumann map, i.e. on the characterization of the unknown 
boundary values in terms of the given initial and boundary conditions. This char- 
acterization employs the same three ingredients introduced in [U] for the analogous 
problem on the half-line: (a) The computation of the large k asymptotics of the eigen- 
functions ^{t,k) and ip{t,k) defining {A{k), B{k)} and {A{k),B{k)} respectively. (6) 
The analysis of the so-called global relation and of the equations obtained from the 
global relation under the transformations which leave invariant the dispersion relation 
of the associated linearized equation, (c) The construction of a perturbative scheme 
for establishing effectiveness. 

This paper is organized as follows: In section [2] we review the main results of 
[8]. In section |3| by employing the ingredients (a)-(c) mentioned earlier, we analyze 
both the Dirichlet and the Neumann problems of the NLS on the finite interval 
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with zero initial conditions. In the former case we express {gi{t) = qx{0,t),hi(t) = 
qx{L, t)} in terms of {go{t), ho{t)}, whereas in the latter case we express {goit), ^o(*)} 
in terms of {gi{t),hi{t)}. Furthermore, in section [s] we show that as L — )• oo, the 
formulas for gi and go respectively, coincide with the analogous formulas obtained 
in [9]. In section |4] we analyze the same problems, but we now express the unknown 
boundary values in terms of the GLM representations. In this connection we correct 
the expressions for {gi(t),hi(t)} in terms of the GLM representations presented in 
[8]. Furthermore, in section |4] we establish the equivalence of the formulas obtained 
directly (i.e. the formulas obtained in section |3]) with the formulas obtained via the 
GLM representations (i.e. the formulas presented in section |4]). 

We emphasize that for problems on the finite interval it is necessary to employ 
all three ingredients mentioned earlier. In particular, the asymptotics of ^{t,k) and 
ip{t, k) yield several possible formulas for the unknown boundary values, so it is ab- 
solutely necessary to employ the ingredient (c) in order to choose the one that yields 
an effective solution. As a warning, we present in appendix [C] a particular set of 
formulas for {gi{t),hi{t)} which do noi yield an effective solution. 

2. Preliminaries 

We consider the NLS equation on the interval [0, L]: 

igt + fe-2A|gp^ = 0, A = ±1, < X < L, < t < r, (2.1) 

where L > is the length of the interval and T > is a fixed finite time. We let (70 (t) 
and hQ{t) denote the Dirichlet boundary values of q{x,t), whereas gi{t) and hi{t) 
denote the Neumann boundary values: 

q{0,t) = go{t), q{L,t) = ho{t), g^(0, t) = 5i(t), qx{L,t) = hi{t). (2.2) 

2.1. Bounded and analytic eigenfunctions. In what follows we present a sum- 
mary of the results obtained in [8]. 



The Lax pair of (2.1) can be written in differential form as 

d (e^('=^+2fc2t)a3^(^^ ^ ^(^^ (2.3) 

where A; G C is the spectral parameter, the closed 1-form W is defined by 

Wix, t, k) = e^{fc-+2'='*)-3 {Qdx + Qdt)fi, as = , (2.4) 



Q=[xq iy Q = 2kQ-i{Qx + \\q\')a:,, 

and a2,A = [a^,A\. We define four eigenfunctions by integrating from the four 

corners of the domain {0 < a; < L, < t < T}: 

fijix,t,k) = 1+ e-'^''''+'^^''^^'Wj{x',t',k), (2.5) 

J{Xj,tj) 



where Wj is the differential form defined in (2.4) with /x replaced by and {{xj,tj)}f 
denote the points (0,T), (0,0), {L,0), and {L,T), respectively, see figure [ij The 
spectral functions {a{k),b{k), A{k), B{k),A{k),B{k)} are defined for /c G C by 

s(k)=(^ ^^^A S(k)=(^^ ^^^^] Sr(k)=(^ ^^^A 
^' \Xb{k) a{k)J' ^' \XB{k) A{k))' ^' \XB{k) A{k)J ' 
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(3) (2^3,^3) (4) 
Figure 1. The contours used to define {fij}f. 



where 

:=M3(0,0,A;), Sik) := fii{0,0,k), S^k) := ii^{L,0,k). 
We also introduce the functions <I>2, Vi) ^2} as follows: 



Let {-Djli denote the four quadrants of the complex fc-plane. Then, 

• a{k) and b{k) are entire functions of k which are bounded in Di U D2. 

• A{k), B{k), A{k), B{k) are entire functions of k which are bounded for k G 

• A;), ^2{t-, k), ipi{t, k), ip2{t, k) are entire functions of k which are bounded 
for k£ D2U Di. 

2.2. The global relation. The eigenfunctions {^j,ipj}f satisfy the following global 
relation: 



c{t, k) = <^2{t, k) b{k)ip2{t, k)e-^^^ * - a{k)^i{t, k)e^^''^ (2.7) 
+ ^i{t,k) \a{k)^2{t,k) - A6(I) V9i (t , k)e'^ikH+2ikL\ ^ ^ 
where c(t, k) is an entire function of k such that 

c{t,k) = Oi—j^ j, /t^oo, feeC. (2.8) 

Indeed, for each t € (0, T), let R{x, t, k) be the solution of the x-part of the Lax pair 
of (2.1 ) such that R{L, t, k) = I, i.e. R is the unique solution of the Volterra integral 
equation 

R{x,t,k) = I + J e^'=(^'-^)^3(Qi?)(x',t,A:)dx, < x < L. (2.9) 

It follows that R is related to /is by 

R{x, t, k) = ^3{x, t, A:)e*'=(^-^)^3^3 ^L, t, k). (2.10) 
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Since 

//3(x, t, k) = fi2{x, t, A:)e-*('=^+2^'*)^^K^), 

this yields 

R{x,t,k) = /i2(x,t,A:)[e-*('^^'+2'='*)^as(A;)]e^^(^-^)^V3H^,i,^)- (2.11) 

Evaluating the (12) entry of ( |2.11 ) at x = and recalling the defin itions (2.6) of 
{<^>j,ipj}l, we find ([27| with c{t,k) := Ri2{t,k). It follows from ([2^ that c(t, fc) is 
an entire function satisfying (2.8). 



In the case of vanishing initial data, the global relation (2.7) reduces to 



(2.12) 



2.3. Asymptotics. Integration by parts in (2.5) shows that 



$i(t,A:) 



k 



+ 



+ 



k 



k ^ oo, k e D2U Da 



<^2it,k) = 1 + 



4^)(t) 



+ 

.(2)/ 



1 



o 



k ^ 00, k e D2LI D4, 



(2.13) 



k'^ 



+ 



-4ik'^t 



— ^ 00, k e D2U D4, 



where 
,(1) 



A; — > 00, k £ D2L) Di, 



9o{t) 
2i ' 

ho{t) 



+ 



+ 



2i 



2i ' ' ' ' ' 4 2i 
and the closed one-form uj is defined by 

A 





$«(t) = 




/ w, 


/ ^, 


Ao,o) 




Ao,o) 


r(L,t) 




r(L,t) 


/ 


/ 


i(L,0) 




i(L,0) 



o;, (2.14) 



a; = - [-ilgpfix + {qqx - qqx)dt] . 



In particular, we find the following expressions for the boundary values: 
go{t) = 2i<P^^\t), 



hoit) = 2t^['\t), 



(2.15a) 



gi (t) = 4$f ) (t) + 2i5o^^'^ (t) , hi (t) = 4^f^ (t) + 2iho^[''^ (t) , 



We will also need the asymptotics of c. 
Lemma 2.1. The global relation {2.1^ implies that 



0<t<T. 

(2.15b) 



c(t, k) 



.(2) 



+ 



+ 



k ' fc2 
^f^(t) , ^S^)(t) + V^S^^(t)(a«+^?)(t)) 



1 



+ 



A:2 



A;3 
+ 



(2.16) 



2ikL 



k^ 00, k € D1UD3, 



where a^^^ = J^^^pjo;. 
Proof. See appendix IbI 



□ 
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Remark 2.2. We emphasize that the proof of lemma 2.1 does not require the knowl- 
edge of the explicit form of c{t, k); it only requires the existence of an entire function 
c{t,k) satisfying ([2^ and 



3. The Dirichlet and Neumann problems 
We will use the following notations: 

• For J = 1, . . . , 4, dDj denotes the boundary of the j'th quadrant Dj, oriented 
so that Dj lies to the left of dDj. 

• dD^ denotes the contour obtained by deforming the contour dD^ so that 
it passes below the zeros of A(A;) in M<o, i.e. below the set { — ^ \ n = 
0, 1, 2, ... }. Moreover, we let dD'^ = -dD^. 

• The functions f+{k) and f-{k) denote the following even and odd combina- 
tions of the function f{k): 

U{k) = f{k) + fi-k), f^{k) = f{k) - fi-k), kGC. 

• A(fe) and are defined by 

A{k) = e^*^-^ - 6-2*'=^, S(A;) = e'^'^^ + e'^'^^. 

The following theorem expresses the spectral functions {A{k), B{k), A{k), B{k)} in 
terms of the given boundary data via the solution of a system of nonlinear integral 
equations. For simplicity, we assume that qo{x) vanishes identically. 

Theorem 3.1. Let T < oo and let qo{x) = for x > 0. For the Dirichlet problem it 
is assumed that the functions goit) and hQ{t), < t < T , have sufficient smoothness 
and are compatible with the initial data. For the Neumann problem it is assumed 
that the functions gi{t) and hi{t), < t < T, have sufficient smoothness and are 
compatible with the initial data. 

Then the spectral functions {A,B,A, B} are given by 



A{k) = $2(T, k) B{k) = -$i(r, k)e'''^ ' , (3.1a) 



A{k) = Lp2{T,k) B{k) = -ipi{T,k)e^''' ^, k £ C, (3.1b) 

where the complex-valued functions {^j{t,k),ipj(t,k)}i satisfy the following system 
of nonlinear integral equations: 

$i(t, k) = f e4^^'(*'-*) [-iAboP^i + (2%o + %)^2] {t\ k)dt', (3.2a) 
Jo 

$2(t, k) = l + \ I [{2kgo - %)$i + i\go?^2\ {f, k)dt', (3.2b) 
Jo 

i-t 

fi{t,k)= / e'^'^^^^'-^'^[-iX\ho\'^(pi + {2kho + ihi)^2]{t',k)dt', (3.2c) 
Jo 

(P2{t, k) = l + X / [{2kho - ihi)fi + i|/io| V2] (*', k)dt' , <t <T, k e C. 
Jo 

(3.2d) 



THE UNIFIED METHOD III 



(a) For the Dirichlet problem, the unknown Neumann boundary values gi{t) and 
hi{t) are given by the following expressions: 



[k^i.{t,k) + igo{t)] 



1 



A{k) 



[kipi^{t,k) + iho{t)] 



+ -^f2-{t,k) + -^—<^2-{t,k) 



+ 



k 



A{k) 



{'P2{t,k) - l)<^>i{t,k)e 



-2ikL 



Mt,k){^2{t,k)-l] 



dk 

(3.3a) 



and 



hi{t) 



T,(k) 1 
^ ^ [k^i-{t,k) + iho{t)] + -—[k'^i^{t,k)+igo{t)] 



JdDO I 2A(A;) 



A(A;) 



+ 



A(fc) 



(^2(t, A:) - l)$i(t, fc) - fe)($2(t, A:) - l)e 



2ifcL 



dk. 
(3.3b) 



(6) For i/ie Neumann problem, the unknown boundary values go{t) and /io(t) are 
given by the following expressions: 



9oit) - 

and 
ho{t) 



1 fm 

7rigooA(/c)\ 2 
+ 



<l>i+{t,k)-^i+{t,k) 



(3.4a) 



$i(t, fe)((/.2(i, fe) - l)e-2^^^ - ($2(t, fe) - l)Mt, k) 
S(fc) 



1 

y^DO A (A;) 



-V9i+(t,A:) + $i+(t,A;) 



(3.4b) 



{Mt, k) - l)«J>i(t, fc) - e2^^^(cl>2(t, fc) - l)</.i(t, A;) 



dk. 



Proof. Equations (3.1) and (3.2) follow from the definitions and (2.5). 



(a) In order to derive (3.3a) we note that the first o f equ ations (2.15b) expresses 

gi in terms of $2^^ and Furthermore, equations (2.13) and Cauchy's theorem 

imply 

-^$(^)(t) = / [<^2- l]dk = f [CD2 - l]dk 

^ J a Do JdDA 



(3.5) 



and 



-^f\t) 



dD2 



Thus, 



i7r$«(t) 



+ 

dD2 JdDi 



[$2 - l\dk 



+ 

dD-j, JdDi 



[k<^l - $f^]dA:. 



(3.6) 



[$2 - l]dk 



(3.7) 



/ [<^2{t,k)-l]dk- [ [^2{t,-k)-l]dk= f ^2-{t,k)dk 

JdDz JdDa JdDz 
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and 



\JdD3 JdDiJ JdDs 



dk 



k<^l - + ^ 



2e- 



-2ikL 



- \ dk + /(t) 



where the function I{t) is defined by 



(3.8) 



m 



r 2e 



(1)1 



dk. 



The last step involves using the global relation (2.12) to compute I{t): 

2Q—'2ikL 



I{t) 



A 



(1) Mil^ ^,Jl)^2ifcL 



k 



+ 



2fc 



- — (/?^^'*] \ dk 



dk 



+ I -^[$i((^2-l)e- 



-2jfcL 



($2-l)v'i]_dA:. 



(3.9) 



The asymptoti cs ( 2.16) of c{t,k) and Cauchy's theorem imply tha t th e first integral 

I /'2^ I I 

on the rhs of (3.9) equals —itt^I {t). Moreover, analogously to (3.7), we have the 
identity 



mip^2\^) — / (p2-{t,k)dk. 

JdD-i 



(3.10) 



Using this identity, the fact that ^p^^ = — 9?2^\ and the expressions (2.14) for ^^""^^ 



and \ we can write the second integral on the rhs of (3.9) as 

I %^2~ - '^[k^i- + ih^] \dk. 
dD° 1 2i A J 



Therefore, equations (3.8) and (3.9) imply 



(go 



A 



[kifi- + iho] >dk 



+ I ^[<I>i((^2-l)e-'*'^^-(«'2-l)H-^^- 



(3.11) 



Using (3.7) and (3.11) in the first of equations (2.15b), we find (3.3a). 



The expression (3.3b) for hi{t) can be derived in a similar way. Indeed, the second 



of equations (2.15b) expresses hi in terms of 1^2^ and (/^^^^ The coefficient Lp^^ is 
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given by (3.10), whereas satisfies the following analog of equation (3.8): 
Jan-, 



kipi - ip\ 



(1) 



2e 



,2ikL 



[k^l-^^l^]\ dk + J{t) 



where the function J{t) is defined by 

[■ ( 2ikL ~\ 



(3,12) 



The last step 

J{t) = 



involves using the global relation (2.12) to compute J{t): 



JdDO 



2 

'a 



kc - <^>\ ' + 



dk 



+ / --[^,i^2-l)-{^2-l)^ie^'''']_dk. 



(3.13) 



asympto tics (|2.16 ) of c(t,k) and Cauchy's theorem imply that the first integral 

le rhs of (3.13) equals —i7Tip\ (t). Moreover, we can write the second integral on 
■hs of (3.13) as 



The 

on the rhs 
the rhs of (3.13) as 



Therefore, equations (3.12) and (3.13) imply 



kT,ipi S (1) 



dk 

(3.14) 



qq in terms 



+ / — [$1(^2-1) -(^2- l)9Jie2^^^] 

Using (3.10) and (3.14) in the second of equations (2.15b), we find (3.3b). 

(6) In order to derive (3.4a) we note that the first of equations (2.15a) expresses 
of Furthermore, equations (2.13) and Cauchy's theorem imply 



-!^$W(t)= / <^^dk= / <^idk 

2 JdD2 JdD4, 



(3.15) 



Thus, 



[ + [ )<^idk= I ^i^dk= I ^<^>i+dk + I{t), (3.16) 



where the function I{t) is defined by 

m = - [ ^ie-''''^^^)+dk. 



10 



J. LENELLS AND A. S. FOKAS 



The last step involves using the global relation (2.12) to compute I{t): 
2 



m 



A 



dD[ 



A 



[$l((^2 



1 e 



-2ikL 



^2Vl\. dk. (3.17) 



The asymptotics (2.16) of c{t,k) an d Cau c hy's t heorem imply that the first t erm on 
the rhs equals —in^^^; equations (3.16), (3.17), and the first of equations (2.15a) 
yield (3.4a). The proof of (3.4b) is similar. □ 



3.1. Effective characterizations. The substitution of the expressions (3.3) for gi{t), 
hi{t) into (3.2) yields a system of quadratically nonlinear integral equations for 
{^j{t, k),ipj{x,t)}i. This nonlinear system provides an effective characterization of 
the spectral functions for the Dirichlet problem. In particular, given the Dirichlet data 
go{t) and ho{t), the system can be solved recursively to all orders in a well-defined 
perturbative scheme. Indeed, substituting into (3.3) the expansions 



i2 



50 = £501 + e^502 H , gi = egii + e^gu + 



ho = ehoi + e ho2-\ , hi = ehu + e hi2 + 

where e > is a small parameter, the terms of 0(e") yield 



(Pj = tpjo + ecpji + e ipj2-\ , j = 1, 2, 

(3.18) 



5in = A / [k^in- + igon] " ^ i^'fin- + iho„] \dk + lower order terms, 

JdD<i 1 2A A J 

(3.19a) 



hln 



S 1 1 

-— [kifin- + ihon] + ^ [k^in- + «50n] >dk + lower order terms. 

(3.19b) 



The terms of 0(e") of the first and third equations in (3.2) yield 

ft 



^lnit,k) 
(pin{t,k) 







^4ife2(t'-t) 



{2kgQn{t') + igin{t'))dt' + lower order terms, 
{2khon{t') + ihin(t'))dt' + lower order terms. 



The odd parts of the latter two equations yield 



$in-(t, k) = Ak / e^*''' 9on(i')c^i' + lower order terms, 

93i„_(t, k) =Ak / e^'''^''^'~^^honit')dt' + lower order terms. 
Jo 



(3.20a) 
(3.20b) 

(3.21a) 
(3.21b) 



It follows from (3.21) that ^i^- and fin- can be determined at each step from the 
known Dirichlet boundary values gon and hon] gin and kin can then be determined 
from ( [3l9l ). 

Similarly, the nonlinear system obtained by substituting the expressions (3.4) for 
go{t), ho{t) into (3.2) provides an effective characterization of the spectral functions 
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for the Neumann problem. Indeed, the terms of 0{e"') of ( 3.4a )-( 3.4b) yield 

90n{t) = - I ^\^<^in+ - ipin+\dk, 

^ JdDl ^ I 2 J 

honit) = -[ ^\ -^ipin+ + ^ln+\dk, 



while the even parts of (3.20) yield 

ft 



(3.22a) 
(3.22b) 



$i„+(t, k) = 2i [ e^*''' c/i„(t')c^i' + lower order terms, 
Jo 

(t, k) = 2i / e'^'-'^^^^'~*^hin{t')dt' + lower order terms. 

JO 

At each step in the perturbative scheme the functions and ipin+ can be de- 

termined from the latter two equations and then g^n and /ion can be found from 



(3.22). 



3.2. The linear limit. The linear limit of (3.3a) yields gi = egu + O(e^) where 

911 = — j I TT^ [k^ii- + igon] - \ [k'fii- + ihon] \dk. 
TT* JdDl 1 2A A J 



Equation (3.21) becomes 

ft 



Thus, 



5-11 



$ii_ = 4A:/' e^^^'(*-*)(7oi(s)rfs, ipu- = ^k [ e^^'='(^-*)/ioi(s)ds. 
Jo Jo 



(3.23) 



This coincides with the formula of appendix [Aj where the linearized equation iqt + 
Qxx = is solved directly. 



3.3. The large L limit. In the limit L — )• 00, the representations for gi and go of 
theorem |3.1| reduce to the corresponding representations on the half-line. Indeed, as 
L — > 00, 

ho 0, hi-)- 0, (pi -^-0, (p2 1, 



A 



— )• 1 as — )• 00 in D3, 



A 



—7- —1 as A; — )• 00 in Di . 



Thus, the L — t- 00 limits of the representations (3.3a) and (3.4a) are 

k 



91 



JdD3 I 2 



-^^--^+'-^<^2-]dk and 50 = -/ ^i+dA;, 

2i 2 J TT Jg£,0 



respectively, and these formulas coincide with the corresponding half-line formulas, 
cf. I9j. 
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4. The GLM approach 



In theorem 3.1 we derived effective representations for {51,^1,50)^0} in terms of 
the eigenfunctions {^j,ipj}f. In what follows we will express the above boundary 
values in terms of the GLM representations. 

For a function f{t, s), we let f{t, k) denote the transform 

ft 



f{t,k) 



The eigenfunctions {<I>j,(/?j}f admit the following GLM representations 

$i(t,A;) = Li-^go(t)M2 + /cMi, ^>2(t, /c) = 1 + L2 + y^oMi + /cMa, (4.1a) 

ifi{t,k) = Ci-UiQ{t)M2 + kMi, ip2{t,k) = l + C2 + Y^oMi + kM2, (4.1b) 

where the functions {Lj{t, s), Mj{t, s), Cj{t, s), M.j{t, s)}i, —t<s<t, satisfy a non- 
linear Goursat system (see [8|) together with the initial conditions 



Li{t,t)='-gi{t), Mi{t,t) = go{t), Ci{t,t) = '-hi{t), Mi{t,t) = ho{t), 
L2{t, -t) = M2{t, -t) = C2{t, -t) = M2{t, -t) = 0. 
Theorem 4.1. Define the function F{t,k) by 

2ikL 



F{t,k) = -ho{t)e''''^M2 - -go{t)M, 



(4.2) 
(4.3) 



+ C2-iX 



hoit) 



Ml + kM2 L 



'-go{t)M2 + kMi 



2ikL 



L2 - iX^Mi + km) Ui - Uio{t)M2 + kMi], 



where C2 is short-hand notation for £2(t, k) etc. Under the assumptions of theorem 
3.1\ the following formulas are valid: 

(a) For the Dirichlet problem, the unknown boundary values gi and hi are given 
by 



9i{t) 



4 



A 



2ik^ 



+ 



k_^ 
A 



Miit,k) 



gojt) 

2ik^ 



k , -2ifcLz7/. L^^ , ^9o{t),.y , %o(i)" 



+ ^(e—F{t,k)).+ 2 



-kM2{t,k) + 



2i 



M2{t,k)}dk, (4.4a) 



hiit) 



2A;2 r,.^ ^ , ^ go(t) 



A;2S r , , , ho(t) 



4 

+ ^iF{t,k)y + '^^kM2it,k) + '^M2{t,k) \dk. 

A 2 2i ' 



(4.4b) 



(6) For the Neumann problem, the unknown boundary values gQ and ho are given 
by 

9oit) = - I XTm {^i^)Li{t, k) - 2£i(t, k) + (e-2^^^F(t, k))+] dk, (4.5a) 
ho{t)=^ [ -^\-^ik)Ciit,k)+2Li{t,k)+F+{t,k)}dk. (4.5b) 
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Remark 4.2. The representations in (4.4) coincide with the representations (4.7)- 
(4.8) in [Sj, except that the last two terms on the rhs of (4.4a), as well as the last two 
terms on the rhs of (4.4b) were missed in 0. These terms, which arise from somewhat 
subtle boundary effects, are needed in order for equations (4.4) to be consistent with 
the representations of theorem 3.1 and are also required in order to obtain the correct 
large L limit. 



Proof. Let us first consider the Dirichlet problem. In view of the GLM representations 

(4.6) 



(4.1), we may write the global relation (2.12 



1 as 



-Li + e'^'^'^ti = kMi - ke^'^^Mi + F - c, 



where 



F{t, k) = -ho{t)e^"''^M2 - ^9o{t)M2 + ((^2 - l)$i - e'''''{^2 - (4.7) 



2ikL I 



The expression of F can be expressed as in (4.3). Letting k — t- —k in (4.6), we find 
-Li + e-'^^^ti = -kMi + ke-^'''^Mi + F{t, -k) - c{t, -k). (4.8) 



Solving (4.6) and (4.8) for Li and £1, we find 
-Li 



2^ 




1 








A 






1 








A 



(e-2^'=^(F-c))- 



(4.9) 
(4.10) 



Multiplying these equations by fee'**^^*^* < t' < t, and integrating along dD^ 
with respect to dk, we obtain 



[ ke^''"^'-''^C,dk 

JdDl 



'^^4ikHt-t')M,dk 

A 



k^^ 



+ I ^e^^^'(-'-''^F.dk, 

JdD<l ^ 



(4.11a) 



A 



k^^ 



A 



e^^^'^'-^'^Midk 



k 



,4ik^{t-t')( -2ikL 



Id DO A 

where we have used that the functions 
k 



{e-'"''^F).dk, 



(4.11b) 



ic{t,k)y 



and 



k 



are bounded and analytic in D^, so that their contributions vanish by Jordan's lemma. 

The next step is to take the limit t' f t in (4.11). This can be achieved by using 
the identities 



ke^'''^^^-^'^f{t,k)dk 



f/(t,2t'-t), 0<t'<t, 



dDi 



f/(i,t), 



0<t' = t, 



and 



k ^iik'^jt-t') 



fit,k)dk 



dDi ^ 



,iik^(T-t') 



f{t,2T-t)dT 



f{t,2t'-t) 
Aik"^ 



dk. 



(4.12) 

(4.13) 
< t' < t. 
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in (4.11 ), we find 



The identity (|4.13|) is also valid if ^ is replaced by k"^ or Utilizing these identities 



71" ^ , , , f 
-Ci{t,2t'-t) = 4 - 

2 JdD'l ^ 



Mi{t,2t' - t) 



dk 



JdD'i ^ Jo 



Mi{t,2t' -t) 



dk 



e^'''^^^-^'^F_{t,k)dk 



and 



vr , , , I k^ 
-Li{t,2t' -t) = A / - 



JdDO A Jq 



Mi{t,2t' -t) 



dDl 



PS 



^Aik^{T-t') 



Mi{t,2T-t)dT 



iik"^ 

Mi(t,2t' - t) 
iik^ 



dk 



dk 



dk. 



Letting t' '[ t in these equations and using the initial conditions (4.2) as well as the 



following lemma, we find the representations in (4.4). 



Lemma 4.3. 



iho{t) 



k 

A 



F^{t,k)dk 



(4.14a) 



+ 



k 



lim / -e^'^^*-^'\e-^^^^F{t,k))^dk 
t'ti JqdI a 



kM2{t,k)dk+ [ ^^^^M2{t,k)dk, 
dDl JdDO 2z 

k 



A 



+ 



kM2{t,k)dk+ , 
dD\ JdDl 2z 



e-^'^^F{t,k))-dk (4.14b) 



kgo{t) 



M2{t,k)dk. 



Proof. We prove (4.14b); the proof of (4.14a) is similar. If one naively takes the limit 



inside the integral in (4.14b), one finds the first term on the rhs of (4.14b). The other 



two terms on the rhs of (4.14b) arise from interchanging the limit and the integration. 



We will next describe how these terms arise in detail. 
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We write 
k 



A 



dDO 



f,^iik-(t-t') I _ ix^Mi ]{Li- ^goM2 ]dk 



+ 



Ml \Mi + M2[ C 



ihr 



k'^M2Mie^"'' dk. 
dDO 



dk 
dk 

(4.15) 



The first integral on the rhs of (4.15) yields the following contribution in the limit 
t' t: 



lim 



igoit) 



t'tt 2 Jg^O 



ke'^''Ht-t')M,(t,k)dk = lim'-^-M2{t,2t'-t) = '-^M^it^t). 

t't* 2 2 4 



On the other hand, utilizing the second row of (4.12), 

igo{t) 



kM2{t,k)dk = ^^^^M2(t,t). 



Therefore, 



lim 



Ml/ ke'^''('-''^M2{t,k)dk='^ 



t'tt 2 Jqj^o 



dDO 



kM2{t,k)dk (4.16) 



+ 



ig^it) 



dDO 



kM2{t,k)dk. 



The first term on the rhs of (4.16) is the contribution obtained by taking the limit 
inside the integral; this term is included in the first term on the rhs of (4.14b). In 
addition to this term, there is also an additional term arising from the interchange of 
the limit and the integration; this is the second term on the rhs of (4.14b). 

We now consider the last integral on the rhs of (4.15), which can be written as 

- / k^~M2Mie^'^''^^'^'^dk = -2 I k^M2{t,k) [\^'''^^^-*'^ Mi{t,2T - t)dTdk. 

JdD° JdDO Jo 



The rhs of (4.17) equals 



-2 / k^M2{t,k)l f e^^^'("-*')Mi(t,2r-t)dr 
JdDO y Jq 



Mi{t,2t' -t) 
Aik'^ 



(4.17) 



dk. (4.18) 



Indeed, the rhs of (4.17) equals this term plus the following expression: 

ft 



-2 [ k^M2it,k)( [ e^''''^^-''h4i{t,2T-t)dT + 
JdDO [J^, 



Mi{t,2t' - t) 
Aik'^ 



dk. 



Integration by parts shows that this integral vanishes by Jordan's lemm a, bec ause 
M2{t,k) is of 0{l/k'^) as A; — ;> oo within D4. Taking the limit t' "[ t in (4.18) and 
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using (4.2), we find that the contribution of the last integral in ( 4.15| ) is 



lim / k^M2Mie^'^^^^-^'^dk 



k'M2Midk + 



aD<l 



2i 



Modk. 



The first term on the rhs is the contribution obtained by taking the limit inside the 
integral. In addition to this term, there is also an additional term arising from the 
interchange of the limit and the integration; this is the third term on the rhs of 



(4.14b) 



Finally, we claim that the limits of the second, third, and fourth integrals on the 



rhs of (4.15) can be computed by simply taking the limit inside the integral, i.e. in 



these cases no additional terms arise. We show this for the term 

„2 



I :-- 



dDO 



the proofs for the other terms are similar. We have 



Co 



We can write this as 

„2 



C2 - —r^Mi 



A 



2 Mij J e^'^^^-'''^Mi{t,2r-t)dTdk. 



f e4^'=^(-*')Mi(t, 2r - t)dr - Mi{t, ^t' -t) 
Jo 4i/c 



dk. 



Indeed, the difference between the preceding two expressions, 



£2 - ^A^i 



A 



e^*'='("-*')Mi(t,2r-t)dr + 



Mi{t,2t' - t) 
Aik'^ 



dk, 



vanishes by integration by parts and Jordan's lemma, since £2 and Aii are of 0{l/k'^) 
as A; —7- c« in D^. Taking the limit t' t, we find 



lim 



JdDO A 



f'^^4ikHt-t')(l^_^_^^^M^dk 



£2 T^Ml 



Mi{t,k) 



gojt) 

2ik^ 



dk. 



However, in this case the additional term 



A 



Co 



vanishes because the integrand is analytic and of 0(1/A;^) as A: —t- 00 in D\. This 



completes the proof of lemma 4.3 



□ 



We now return to theorem 4.1 and consider the Neumann problem. Solving (4.6) 



and (4.8) for M\ and A^i, we find 



kMx 





S 1 















(4.19) 
(4.20) 
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Multiplying these equations by /ce^*'^^^* < t' < t, and integrating along dD\ 
with respect to dk, we find 



-Mi{t,2t' -t) 



dDl 



vr 



Mi{t,2t' - t) 



dDl 

where we used that the functions 

are bounded and analytic in the interior of dD^ so that their contributions vanish 
by Jordan's lemma. Letting t' t m. these equations and using the initial conditions 
(4.2) as well as the following lemma, we find the representations in (4.5). 

1 



A 



A 

2 - 1 



-2ikL 



A 



A 



A 



F)+ }dk 



{e-^"''^cit,k))+ 



Lemma 4.4. 

lim [ ^e^''''^'-''^F+(t,k)dk 



A 



F+{t,k)dk, 



(4.21a) 



lim 



1 



t'tt Jqdo a 
Proof. Note that 



^''''^'-''\e-^'''^F{t,k))+dk 



1 



{e~^'''^F{t,k))+dk. (4.21b) 



J dDl A J dDl A [ 2 



S(L2-a^Mi) {Ci-^ho{t)M2 



- AkM2 ( A - -ho{t)M2 J - A (^La - iAy Mi) fcTWi - Sfc^Mz^ii [> . 



Equation (4.21a) now follows by using arguments similar to those that led to lemma 



4.3 The proof of (4.21b) is similar. This completes the proof of lemma 4.4 and hence 



of theorem 14.1 



□ 



4.1. Equivalence of the two representations. We will show that the representa- 
tions derived using the GLM approach in theorem 4.1 coincide with those of theorem 



4.1.1. The representations for gi and hi. Using the expression (4.7) for F as well as 
the formulas 

-^i = ^V'i-, j = l,2, 



^ 2A; ^ ' 



we can write the representation (4.4a) of 51 as 
4 



9i{t) 



J dDl 



^ [k^i^ + i/io(t)] + ^ [k^i-{k) + i9o{t)] 



+ 



k 
A 
k 
A 



(992 - l)^ie 



'2ikL 



{^2 - 1)^1 



4k 



, ^9o{t)^ , go{t) _ , 
+ ^^2- + -^V2- }dk, 
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The identities 



JdDO 



ip2~{t, k)dk = —7riip^2^ = tti^p^ 



.(1) 



<f2-dk, 



imply that the function ip2- in the above integrand can be replaced with ip2-- More- 
over, since the integrand is an odd function of k, the contour dD^ can be replaced 
with dD^ In view of the identity 



k 

'a 



4k 



igoit) 



$2 



we find the representation for gi in (3.3a). Similar computations show that the 
representations for hi are also equivalent. 



4.1.2. The representations for and /iq. Using the expression (4.7) for F as well as 
the formulas 

$1+ = 2Li - igoM2, 'Pi+ = 2Ci - ihoM2, 
a straightforward computation shows that the representations (3.4) and (4.5) are 
equivalent. 

Appendix A. The linear limit 

In this appendix we analyze the linearized version of the NLS equation, iqt+qxx = 0, 
on the interval [0,L] using the unified method of The global relation for this 
equation is (cf. equation (2.8) in [7j) 



qoik) - ~g{ik^) + e-'^^h{ik') = e^'^' grik), k e C, 



(A.l) 



where 



Qo{k) 



-ikx 



qo{x)dx, 



qrik) 







-ikx 



q{x, T)dx, 







g{k) = -kgo{ik^) + igi{ik^), 

rT 

'„2\ _ I Ak'^s 



9j{ik 



h{k) = -khoiik'^) + ihi{ik'^), 
rT 



gj{s)ds 



hjiik"^) 



„ik s 



hj{s)ds, j = 0, 1, 



and gj{t),hj{t), j = 0,1, denote the Dirichlet and Neumann boundary values as in 
(2.2). Equation ( |A.l ) and the equation obtained by letting /c — )• — A; in ( A.l ) are the 
following equations: 



qoik) + kgoiik"^) - m{ik^) " kho{ik^)e-'^^ + ihi{ik'^)e-'^^ = e'^^^ qrik), 



- kgo{ik^) - igi{ik^) + khQ{ik'^)e^^^ + ihi{ik^)e 



„2\ikL 



Jk^T 
Jk^T 



qri-k). 



Eliminating hi from these equations and then solving for gi, we find 



giiik^) 



A{k/2) L 



e^'^'^ie'^^qTik)). 



AkL 



qo{k))- - kJ:{k/2)go{ik'^) + 2kho{ik'^ 



Multiplying this equation by ke "^^^^ and integrating along dD\ with respect to d/c, 
the term involving qx is eliminated and we find 



ike 



{e'^^qoik))^ - k^{k/2)go{ie) + 2kho{ik' 



Iqdo A{k/2) L 

Assuming that qo = and performing the change of variables k = —21, we find 



9i{t) 



2i^j:{i) 



go{s)ds- l\^''"^'-'ho{s)ds 



dl. 
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Using the identity 

fc2 rT 



dDO A{k) 



Jo 



K{t,T) 



(A.2) 
dk, 0<t<T, 



and a similar identity obtained by replacing by in (|A.2|), we arrive at 



9iit) 



4 f r2/2s(/) 



i A(/) io ' 
Ail) 



gois)ds 



m)9o 
2iA{l) 



dL 



which coincides with formula (3.23) for gn. 



Appendix B. The asymptotics of c{t,k) 



We will prove lemma 2.1 We will make no assmnption on the function c{t, k) except 
that it sa tisfie s the global relation ( |2.7[) and that it has the boundedness properties 
stated in (2.8). The functions {^j, satisfy the systems 



and 



^it 

^2t 



'Pit 

V2t 



-4i/c2$i - iX\go\^^i + (2%o + %)^2, 
A(2/c5o -%)^i + Ai|5oP^2, 



— 4iA;^(/3i — iA|/iop9?i + {2kho + ihi)ip2, 
X{2kho — ihijipi + Ai|/ioP</'2- 



In view of the initial conditions 

$i(0,fc) = (/Pi(0,A;) = 0, 
this leads to the asymptotic expansions 



$2(0, A:) = (^2(0, A:) = 1, 



^i{t,k) 
^2{t,k) 



+ 



f ^iit,k) 
\^2{t,k) 



+ 



1 



k 



+ 



^^^\t) k^ 



+ 



cI>W(0) 



+ 



(o) + ^S'^(o)/i:o> 
-fgom^\o) 



(0,t) 



1 

P 
1 



+ 



+ 



+ 



A;2 +^(a;3 



.{L,t) ^ 



1 



o 



(B.l) 



(B.2) 



k — )• 00, k G 



00, k £ 
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Similarly, we have 

fb{k)\ fO 



\a{k) 
+ 



+ 



5(i)\ 1 /6(2)\ 1 



k 



+ 



(B.3) 



2ikL 



k — )• oo, k £ C 



Substituting these expansions into the global relation (2.7), we find 

1 



c{t,k) 



O 



+ 



fc2 



2ikL I Aik^t 



+ 



6(1) _$W(0) 



k 



+ 



k^ 



+ 



qo{L) 
2i 



+ 



2ikL 



ci>«(t) ^ c|>f)(t) + .|>W(tKaW + 4^)(t)) ^ ^ 



k ' k"^ 



(B.4) 



k 



+ 



A:2 



+ 



k — >• cxD, A; G 



Th e ass umption that c{t,k) is of 0((1 + e )/k) as A; — )■ oo implies that the terms 
in (B.4) involving e 



-iik^t 



and e^**^^* must vanish, i.e. for consistency we require 



6« = <I>S^)(0), 



*fl = v.S^)(o). 



(B.5) 



Using the expressions for b^^\ <I>^^^(0), and ip^^\o), we find that equations (B.5) are 
valid iff the initial and boundary conditions are compatible, i.e. iff 



qo{0) = 5o(0), qo{L) = ho{0). 



Appendix C. A non-effective characterization 



The representations of gi and hi derived in theorem |3.1| yield an effective char- 
acterization of the spectral functions {A{k), B{k),A{k),B{k)} in the sense that the 
resulting nonlinear system can be solved uniquely at each step in a well-defined per- 
turbative scheme. 

In what follows we present a set of non-effective formulas (cf. Eqs. (14)-(16) in 
[3]). Let 



R{0,t,k) 



d{t,k) c{t,k) 
Xc{t,k) d{t,k) 



where R is defined in (2.9). It follows from (2.10) that R satisfies 

Rt = -2ik^[a3,R] + Q{0,t,k)R- Re'''^^''Q{L,t,k). 
This implies that {c, d} satisfy the system 

,2 



ct = —Aik c ■ 



Ai(|5o|^ + \ho\^)c + (2%o + i9i)d - {2kho + ihi)e^^^^d, 



dt = Ai(|go| ~ l^ol )d + X{2kgQ — igi)c — \{2khQ + ihi)e 
c{0,k) = b{k), d{0,k) = a{k), 



2ikL- 



(C.l) 



(C.2a) 
(C.2b) 
(C.2c) 



THE UNIFIED METHOD III 



21 



where c and d are short-hand notations for c{t,k) and d(t, k ), respectively. Let us 
consider the Neumann problem. Integration by parts in (2.9) shows thalj^ 



c{t,k) 



|^ + 0(A;-2) + 0(e2*'=^fc-i), 



k G C+, 
k G C-, 



oo. 



This gives the representations 

1 
vr 
1 
vr 



goit) 

ho{t) 



c+(i, k)dk, 

[cos{2kL)cj^{t, k) — i sin(2A;L)c- (t, k)]dk, 



(C.3a) 
(C.3b) 



which can be used to eliminate {^Oi^o} from (C.2). The resulting system for {c,d} 
is formulated only in terms of {gi, hi}. However, this system is not effective. Indeed, 
substituting into (C.2)-(C.3) the expansions 

i2 



d 



1 + e^do + 



eci + e C2 H , 



the terms of 0(e") yield 



e^""^^ \2kgon + igin - (2A;/io„ + ihin)e'^'^^] + lower order terms (C.4) 



and 



90n 



^n-^dk^ 



[cos(2A;L)c„+ — z sin(2/cL)c„_]dA;. 



Thus, in order to determine {gon, hon} at each step of the perturbative scheme, we 



need to know the combinations Cn+ and cos(2A;L)cn+ — isin(2A;L)cn-. Equation (C.4) 
shows that these combinations satisfy 

t 

+ lower order terms. (C.5b) 



'^e"" Cn+}^ = e^*^ [2igin — 2khon^{k) — ihin^{k)] + lower order terms, (C.5a) 
{e^^^' [cos(2A:L)c„+ - i sin(2A;L)c„_] }^ = e^^^' [-2kgonA{k) + %„S(A;) - 2«/ii„] 



Although the function gon has been eliminated from the rhs of (C.5a), the function 
hon remains unknown. Similarly, although the function /iQn has been eliminated from 
the rhs of (C.5b), the function gQn remains unknown. This shows that the solution is 
not effective. 
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